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ABSTRACT 


We have used computer models of the microwave emission 
from the earth ' s atmosphere to study the problem of retrieving 
meterological information from the SMMR instrument that will 
be flown on NIMBUS-G. We describe methods of retrieving rain 
rate, wind speed, cloud height, and ocean temperature when the 
satellite is over the ocean. 
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FOREWORD 


This report describes the work performed by Systems and 
Applied Sciences Corporation (SAS) under NASA Contract 
No. NAS5-24191. The work involved analysis of the problem of 
determining certain metero logical parameters from the Scanning 
Multichannel Microwave Radiometer (SMMR) instrument aboard 
NIMBUS-G, when the satellite is over the ocean. To perform 
this analysis, we used a computer code developed by SAS under 
a previous contract to generate models of the microwave emission 
from the earth's atmosphere. The models specifically include 
the effects of rain drops. 

An ensemble of models was generated that simulate the 
microwave emission under various conditions of rain rate, ocean 
temperature, cloud height, and wind speed. Different methods 
were used to retrieve these parameters from the calculated bright- 
ness temperatures. We evaluated the accuracy, of each method, 
and determined the limitations on our abilities to retrieve the 
above mentioned meterological parameters from the SMMR data. 

The balance of this report describes in detail the methods 
used in the retrievals and detailed estimates of the accuracy 
attainable with the SMT4R instrument. 
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REMOTE SENSING OF RAIN OVER THE OCEAN 


SECTION 1 - INTRODUCTION 


Data from the Scanning Multichannel Microwave Radiometer 
CSMMR) aboard NIMBUS-G can be used, when the satellite is over 
the ocean, to determine wind speed, rain rate, cloud height, 
and ocean temperature. The experiment consists of a ten-channel 
radiometer that observes both vertical and horizontal polariza- 
tions at 6.6, 10.69, 18, 21, and 36.9 GHz This report describes 
a theoretical investigation into the retrieval of the four above- 
mentioned meterologial parameters from SMMR observations. The 
investigation involved constructing numerical models of the 
microwave emission from the atmosphere above the ocean and in 
the presence of rain and varying the parameters of the models to 
generate an ensemble. The ensemble is then used to generate 
regression coefficients that give the value of the parameters in 
terms of the brightness temperatures. By doing this, we can 
investigate the probable limits of our ability to retrieve the 
desired parameters, and devise the best strategies for doing so. 

We adopt, for our model, an atmosphere that is horizontally 
homogeneous. The lower boundary is the ocean surface. The micro- 
wave emission arises both from the atmosphere and the surface of 
the ocean; in addition, the ocean reflects microwaves. Because 
the reflection coefficients are different, in general, for the 
horizontal and vertical polarizations, the radiation reflected 
from the ocean is polarijed. The effect of rain on the microwave" 
radiance is both to increase the intensity at every wave- ! 

length and reduce the degree of polarization. The latter effect 
comes about because, as microwaves- are scattered by the raindrops, 
the plane of polarization changes. 
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SECTION 2 - MODELS 


The computer models that were generated for this 
investigation have four input parameters: rain rate, cloud 
height, surface wind speed, and ocean temperature. Other para- 
meters that may in -fact vary from place to place (or time to 
time) have been fixed, to make the problem tractable. 

These models are used to simulate the observations under the 
conditions where the meterological conditions do not vary within 
the field of view of the SMMR instrument. Additional models 
were made where only part of the field of view is covered with 
rain; this is discussed in Section 5. 


The models have a constant vertical temperature lapse rate 
of 6.5 K .km with the 273.1 K level fixed at the cloud top. 
The relative humidity varys linearly between 80% at the ocean 
surface and 100% at the cloud top. We have a layer 0.5 km 
thick, with its upper boundary at the cloud top, that contains 

-3 

non-precipitating liquid water with a density of 0.5 g ra . 


The sources of opacity are molecular oxygen^ water, vapor, 
the non-precipitating liquid water, and the raindrops. 

We consider the size distribution of the raindrops to be 
given by the Marshall-Palmer distribution 

,N(r) dr = .16 e ^^'dr cm ^ 

- 0.21 

where r is the drop radius (cm), and A = 81.56 R ' , where 

R is the nominal rain rate in mm hr ^ (Wilheit et al., 1976, 
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and references cited therein) . The scattering cross-section 
is calculated for each drop size separately according to Mie 
theory and the composite scattering function is found from 


S{0,0' ) = X)N(r) p(0,0' ,r) 

where S(0,0') is the total cross-section for scattering from 
an angle 0 to an angle 0', and p(0,0',r) the cross-section 
for a sphere of radius r. 


The reflection coefficient of the ocean surface is a 
function of temperature, wind speed, and frequency. The re- 
flection coefficient is calculated according to a model by 
Wilheit (1978) . For wind speeds less than 7 m s the 
reflection coefficient is linear in wind speed; it is also 
linear for wind speeds greater than 7 m s There is a 

discontinuity in the slope at 7 m s ^ that is due to the for- 
mation of surface foam at that wind speed. 


The solution of the radiative transfer equation follows 
the method of Herman (1963). The atmosphere is divided into 
a number of layers with equal optical depth, and the equation 
of radiative transfer is integrated in each layer for each of 
12 polar angles. The solution is done by iteration; fewer than 
10 iterations are required for convergence. We note that, 
because there is no incident radiation field at the top of the 
atmosphere (the 3 K cosmic background is negligible) , the radia- 
tion field is azimutually symmetric. This symmetry reduces the 
computing time by about an order of magnitude over what would 
be required in the more general case. 
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Two aspects of the models that are fixed are the relative 
humidity profile and the amount of non-precipitating liquid 
water below the freezing level. At each frequency, the optical 
depth of the atmosphere excluding the contribution of rain- 
drops is fixed by the choice of freezing level. For reference 
purposes, we have included the Appendix A table showing thei 
total water vapor (cm) and the non-rain optical depth at each 
frequency. 
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SECTION 3 - RETRIEVAL 


We generated an ensemble of 1880 models for use in the 
retrieval studies. In addition to the four independent para- 
meters in the model — rain rate, freezing level, wind speed, 
and ocean temperature — it was found advantageous to add a 
fifth parameter, the scattering optical depth at 6.6 GHz, to 
the set of parameters to be retrieved. This last parameter, 
of course, is not independent of rain rate and freezing level. 


The retrievals in this section deal only with what we 
call the "homogeneous" case, where the meterological parameters 
are constant within the field of view. The "non-homogeneous" 
case is treated in Section 5, after we have discussed the 
general retrieval method and the accuracy of.. the (somewhat . 
artifical) homogeneous retrievals in Section 4. 


To do the retrieval, we represent each meterological 
parameter p^ as a linear combination of functions f^ of the 
10 brightness temperatures T ^ . In practice, we used 3 sets 
of functions 


Linear 


f. = T. 
3 3 


Quadratic 


f. = T. 

3 3 

f.+lO = T.' 
3 3 


Logarithmic f ^ = £n 


T . 
3 


j = 1,10 


where the a. are some constants chosen so that a. - T. > 1. 

1 3 3 
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We wish to represent each parameter as 




P . - ^ a . . f . 

1 J J 


and determine the a. . so that the residual be minimized in the 

13 

least-square sense. We must consider two sources of error: one 

due to the difference between the fitted and the true functions, 

the other due to noise in the observed brightness temperatures. 

2 th. 

Let Oj be the noxse in the j — channel; the error in the retrieved 

value of p. due to receiver noise will be 
^1 


Ni 



3 


a . 2 

13 


9f. 


if f . are chosen so that the average valve of : 
3 


is 1. 


The error in the fitting function will be 


2_1 
’^i N 



and we seek to minimize the function. 


2 2 

+ Si 


~ Ea. • 

"11-1 N " ki ^ XI 


<2 1 


ij ^kj‘ 


th 

Here, p . is the value of the i— parameter in the k model, 

t h 

f, . the value of f . in the k — model, and N is the number of models 
k3 3 
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We determine the covariance matrix S . . by 



where the superior bar represents the average over the 

th 

we define the covariance between the i — parameter and 
by 


ensemble ; 

f . by 
D 



We now take the derivitive of X^+ within respect to each 

and set it equal to zero; this system of equations gives the 

values of . that minimize X? + • The solution is 

13 Ni 


a 


ii 5 *^ik ^kj 


-1 


where ^ ^ is the inverse of the matrix 


S. . = S. . + D . . 

ID ID 3D 


. . th • 2 

where is a diagonal matrix whose jj — term is a ^ . 


This system can be solved either by matrix inversion or by finding 
the eigenvalues and eigenvectors of ^ , and solving for the ^ in 
this eigenvector basis. We have chosen the latter approach, since 
the eigenvalues of give us significant information about the 

nature of the retrieval problem. Consider, for the moment, the 
case where all the a^. are equal and uncorrellated. 
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In this case, is a scalar multiple of the identity matrix; 

y\ 

S.. =S.. +0^ 6... 

ID 13 ID 

y\ , 

The eigenvectors of are the same as the eigenvectors of , 

and the eigenvalues A, of S. . are related to the eigenvalues 
/N ^ k ij 

^ of by 

Let be the k— eigenvector of (or ; we v/ish to represent 

(for the moment) 


Pi 


a' . . 6 . 
y ID D 


where 3 • 
D 




n 


■jm m 


Then a ' . , = 

Ik 


? Y 0. . e 

D 3 -D Dk 


and it is clear that if | | is very small 

and the contribution to 



will be very large 


2 
Wi 




will be correspondingly very large. 


We must, therefore, keep in mind the values of 0. when we design 

2 ^ 

the retrieval, for tv 70 reasons. One is that, when 0 = q.., we would 
have 

L i 

j 


a 


ik 


A. 


0 . . e 
ID Dk 


and the expansion coefficients would tend to infinity for small 
values of A . The other reason is that eigenvectors with values of 
Aj^ < 0^ contribute very little information to the retrieval. We 
wish to retrieve the values of four parameters; we can do this, 
roughly speaking, only if at least four of the A^^ are greater then 
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We should, at this point, bring up the question of the 
validity of a linear retrieval method, since the brightness 
temperatures are not linear in the parameters. We shall 
reserve judgment on the utility of non-linear methods until 
we have discussed the following point. 

The field of view of the SMMR instrument, at the longest 

wavelength, is 156 km on a side. We can be certa.in^' .a priori , 
that the meterological parameters will not be constant within the 
field of view. We wish to retrieve average values of the para- 
laeters, and want these values to be unbiased. In this context, 
unbiased means that the average of a large number of values of 
some parameter p retrieved from the observations will not be 
systematically different from the true average value of that 
parameter. 


Let us break up the field of view into N regions of equal 
area, with the conditions uniform in each region. We take, as 
an example, the case where we retrieve the parameters as 


Pi 


E 

j 


la . . T . 
11 1 


. B.. 




• th 


where is the average value of the j 

over the field of view. Let T, . be the brightness temperature 

. , . th , , . , , , th ^ . 

in the J — channel in the k — region, 

+- v» 

i— parameter in that region. The value we wish to retrieve 
is 


brightness temperature 

ness temperature 

and p, . the value of the 
ki 


<p.> = 


N 


E 

k 


ki' 
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However, the value we will retrieve is 



The difference between <P.> and P., which is a bias in our 

1 1 

retrieval, is 


p.-<p.>=-^ 3-- 

^ j 


E -E 


kj 


we note that the term in square brackets is the variance 
if Tj over the field of view. Thus, we can write 


<Pi> = - E . var 

j ■' 



We see that by including quadratic terms in the retrieval, we 
may increase the accuracy of the retrieval for the case where 
the conditions are homogeneous within the field of view at the 
expense of biasing the result when the conditions are not 
homogeneous. This point will be treated again in Section 5. 
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SECTION 4 - ACCURACY OF THE RETRIEVALS 

In the section, we shall consider the (somewhat unrealistic) 
case where the meterological parameters are constant within a 
field of view of the SMI4R instrument, which we call the "homo- 
geneous” case. The problem introduced when the rain rate is 
allowed to vary within the field of view (FOV) ( "non-homogeneous" 
case) will be treated in the next section. The primary purpose 
for presenting the results in this section is not one of 
developing a retrieval algorithm — for this, we must consider 
the non-homogeneous case — but to provide a comparison of the 
non-homogeneous case, and underscore the importance of consider- 
ing the case where the rain rate varies over the FOV. 

We have written a program that used the ensemble of 
calculated brightness temperatures to calculate expansion co- 
efficients for retrieving freezing level, rain rate, wind speed, 
sea temperature, and the optical depth at .22 cm The pro- 
gram then calculates the standard deviation of the retrieval, 
including the effect of the noise in the radiometer. In Table 1 
we present some representative results for the homogeneous case, 
where we have made a linear, quadratic, or logarithmic fit to 
the ensemble of models. Columns 1 and 2 show the range of x 
over which the subset of the entire ensemble was selected. 

Column 3 shows which type of retrieval was used. Column 4 shows 
the restriction on the v/ind speed in the models selected, if any. 
The next five columns give the standard deviations in freezing 
level (h) , rain rate (r) , wind speed (v?) , sea surface temperature 
(T^) , and rain opacity (x). For the quadratic retrievals, the 
biases discussed in Section 3 are also given. 
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We see that, as we v/ould expect, the error in the retrieval 
increases as t (and hence the rain rate) increases. We can also 
see that the quadratic retrieval is more accurate than the 
linear retrieval because it fits the relationship between the 
brightness temperatures and metero logical parameters better. 

The logarithmic retrieval is no better than the quadratic 
retrieval, is harder to analyze for bias, and will not be con- 
sidered further. 

In Table 1 we show the bias in a quadratic retrieval, as 
we discussed in the last section- While we cannot know what the 
variance in T„ will be over the field of view, for computational 
purposes we have used the variance of each T^ in the subset of 

IJ 

the ensemble that was used to calculate the retrieval coefficients 
The biases calculated in this manner are not small compared with 
the standard deviation of the retrieval; this supports the notion 
that we should not use quadratic terms in the final evaluation 
of the parameters, but rather use a linear fit over a suitably 
restricted range of t. VJhich range to use might be decided by 
using a quadratic fit over the entire range of t to find a 
provisional value of t . 

We also note from Table 1 that for some of the retrievals, 
the subset of models used to calculate retrieval coefficients 
was restricted to wind speed greater than 7 m s ^.or to wind 
speeds less than 7 m s ^ (the relationship between speed and 
ocean reflectivity has a discontinuity in its slope at 7 m s ^) . 

We see that, by using coefficients taken from the appropriate 
subset of models, and having an approximate value for wind speed, 
we can make a modest improvement in the retrieved value of the 
wind speed, and a considerable improvement in the value of the 
ocean temperature - 


4-2 



The emissivity. of. sea. water, changes due to changes, -in both wind 
speed and sea surface temperature. Having some a' priori 
knowledge of wind speed helps to. reduce the uncertainty in the 
sea temperature. 

We now turn to a different problem. Because only one 
antenna is used for all five frequencies, the field of view is 
different at different frequencies. Thus, at the lowest fre- 
quency, the effective integration time is longer than for higher 
frequencies, and therefore, the receiver noise will be smaller. 
This fact somewhat complicates the analysis of the retrieval 
problem. First, let us look at the effects of deleting some 

channels from the retrieval, but leaving the receiver noise 

2 

the same for all channels (here a = 0.04 K) . In Table 2 we 
show how the standard deviation in each parameter changes as 
we delete different channels. We see that, with only a small 
increase in the error, the wind speed and ocean temperature can 
be found from the 6.6- and 10- GHz channels alone. This cor- 
responds to finding the average values over the largest field of 
view for these parameters; using the higher- frequency channels 
will, not give significant information about the wind speed and 
sea temperature on a smaller spatial scale. 

For the data analysis, four cell sizes have been adopted 
that correspond roughly to the field of view at different 
frequencies. The cell sizes are 156 km at 6.6 GHz, 97.5 km 
at 10 GHz, 60 km at 18 and 21 GHz, and 30 km at 37 GHz. We can 
retrieve the wind speed and ocean temperature only on the 156 km 
scale. 


We have a choice of how to proceed with the freezing 
level and rain rate. These parameters can be retrieved from the 
18- , 21- , and 37- GHz data alone, on a 60-km scale. It probably 
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would not help the retrieval much to include the 6.6- and 10-GHz 
data from the appropriate 160-km cell, as, while the retrieval - 
is somewhat improved by adding the low-frequency channels, the 
average value over the 160-km cell may not be representative of 
the conditions in a particular 60-km cell. Similarly, we could 
try to calculate rain rate and freezing level on a 30-km scale, 
but this would require using average values on a 60-km cell for 
the 18- and 21-GHz channels. It is not clear without further 
study and some experience with real data what the retrievals on 
the 30-km scale would mean, although analysis of the 37 GHz data 
on this scale certainly will give information on how uniform 
conditions are on the larger scale. 

In Table 3 we show how the standard deviation in the 

retrieval changes as we use different cell sizes. The values 

2 . . 
of a used are gxven on the left, and the standard deviations of 

the retrieval of the parameters are given on the right. The 
reason that the noise changes with cell size is that the inte- 
gration time is longer if, for example, we average the 37 GHz 
data over four 30-km cells than it is for one 30-km cell. The 
noise figures adopted here are believed to be representative of 
the SMMR instrument (which has not been launched at the time of 
this writing) . In each group, there are four cases presented, 
that correspond to retrievals on the 156- , 90- , 60- , and 30-km 
cell sizes. For the 90- and 60-km cells, we delete the 6.6-GHz, 
and the 6.6- and 10- GHz channels respectively. The fourth re- 
trieval, on the 30-km scale, uses the 37- GHz. channel for the 
30 km cell, and the 18- and 21- GHz channels from the correspond 
ing 60- km cell. 

2 

In our calculations, we have calculated the values of X- 

2 . ^ 
and a separately for each parameter. In the case of a linear 

retrieval over some range in t , is reduced by reducing the 
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2 

the size of the interval, while is not appreciably affected 

by the size of the interval. In an attempt to develop an optimal 

retrieval algorithm, it would seem that we should divide the 

2 2 

range of x into intervals small enough so that ~ and 

then compute a different set of linear retrieval coefficients 

for each interval of x. Then, for each observed set of ten 

values of T^, we would decide which interval the value of x is 

in and use the appropriate set of retrieval coefficients. In 

this sense, the size of the interval will depend on the values 

of o . . 

1 


However, as we shall see in the next section, there are 

other considerations that will preclude making the intervals 

2 2 

in X small enough so that y, - o.. It seems, for our 

2 

analysis of the retrieval problem, that X- will be larger than 

2 ^ 

0^, and as long as the noise values that obtain in the space- 
craft are not too much larger than the ones v;e have used, the 
results will not be significantly dependent on these values. 

The largest errors, by far, in the retrieval are due to the 
non-linear nature of the dependence of the brightness temperatures 
on the meterological parameters and to the uncertainties of 
calibrating the spacecraft observations with suitable ground 
truth. 
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SECTION 5 - PROBLEMS WHEN THE FIELD. OF VIEW 
IS NON-HOMOGENEOUS 


So f:r, everything we have said concerns only the case 
where the meterological conditions are homogeneous within 
the field of view (FOV) . We must be prepared, however, to 
deal with the situation — which probably will occur more 
frequently than not — where only som.e fraction f of the FOV is 
covered by rain. To investigate this problem, we constructed a 
set of what we call "non-homogeneous" models. We use the homo- 
geneous ensemble to generate the non-homogeneous set by, for a 
given sea temperature, cloud height, and wind speed, constructing 
a set with 

^Bi = f ^Bi(^) + (1 - f) T3.IO) 

r'= fr 
t'= ft 

til 

where T . (r) is the brightness temperature in the i channel 

-1 

at the rain rate r, x is the raindrop opacity at .22 cm , 
and 0 ^ 1 . 


We can calculate the retrieval coefficients from either 
the homogeneous or non-homogeneous ensemble of models. For the 
latter case, we also calculate the coefficients to retrieve the 
fractional area covered by rain, f. (In this case, we do not 
calculate coefficients to retrieve the cloud height.) We call 
these sets of coefficients, again, "homogeneous" or " non-homogeneous " . 
We have used both sets of coefficients to retreive the parameters 
of both ensembles of models over several ranges of x. Some 
results are shown in Table 4 . The first two columns show which 
coefficients and which ensemble was used (H = homogeneous, 

N = non-homogeneous) . The subset used for the particular 
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retrieval was based on the value of t ; the smallest and largest 
rain rates and the number of models are given in the last three 
columns. Tables 5 and 6 compare retrievals based on all 10 channels, 
and only the 6 highest-f requency channels, respectively. The 
shov7 the standard deviation in the retrievals; f is the fraction 
of the FOV covered with rain, and its retrieval is only valid 
when both the ensemble and coefficients are non-homogeneous . 

We see, as we did in Table 2, that omitting the 6- and 10- GHz 
channels has only a small effect on the retrieval of f and the 
rain rate r. 


The first point to note is that we can, with some confidence, 
retrieve the value of f if we have used the non-homogeneous 
coefficients with the non-homogeneous ensemble of models. The 
standard deviation of the retrieval of f is about .02 - .08, 
while 0 < f < 1 . 


The second — and more important — point is that we cannot 
use the homogeneous coefficients with the non-homogeneous 
ensemble of models. For example, if we use the H coefficients 
with the N ensemble of models, in the range .01 < t < .05 
(2.25 < r < 15) we find that standard deviation for rain rate 
is 9.1 mm hr ^ for the linear retrieval, which is a useless 
result. 

It is also clear that retrievals using quadratic terms 

fare very poorly; the reason is that the importance of the 

2 2 2 . ... 

difference between (T^^ + T 2 ) and T^ + T 2 is magnified in 
the non-homogeneous case (see the discussion of the bias on 
§III). 
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This illustrates one difficulty of relying on models to 
interpret the data from the SMMR instrument — or any other 
sounding instrument. V7e have just seen a retrieval scheme,' 
as we outlined in Section 3, that looked good with computer 
simulations, which would in fact have been very bad- We also 
must note, here, that the problems will increase even further 
if, i) there are more than two rates of rain in the field of 
view, and; ii) the other parameters vary also. 
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SECTION 6 - COMJIENTS 


A few words should be said about the problem of verification 
of the SMMR retrieval. In view of what has been said above, it 
is clear that the SMMR data will be worthless without verification 
of the retrieval values of the meterological parameters by 
comparison with ground-based observations. For the moment, let 
us consider only the rain rate problem. 


Common experience tells us a few things about the nature 
of the rainfall distribution. One is that the spatial rainfall 
pattern is more or less constant over a time scale of about 
one hour. The second is that, as we view it from one point on 
the ground, the rainfall rate usually changes very rapidly on 
this time scale. This rapid change is due, in large part, to 
the motion of the rainfall pattern over the ground. Thus, the 
common experience is that, on a scale of 30 - 150 km, the rain- 
fall rate will not be uniform most of the time. This qualitative 
argument alone suffices to indicate that we must consider the 
retrieval problem in the case where only some fraction f of the 
field of view is covered by rain clouds. 


One type of observation potentially available to us is data 
from the aircraft SliMR simulator. Since the aircraft altitude is 
very small compared to that of the satellite, the field of view 
is corresponding smaller. Hence, experiments done with the air- 
craft underflying the satellite, and sampling rapidly many points 
in the satellite field of view, v/ould give us information as to 
the distribution of brightness temperature in the satellite 
field of view. 
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The problem, of course, is that the speed of the sub-satellite 
point for NIMBUS G will be about 2300 km hr which is somewhat 
faster than any aircraft can fly. Even getting a minimal sample 
of aircraft data would require, probably, more flight time than 
is feasible. While it is probably not feasible to use aircraft 
data to calibrate the SMMR data, such observations would be 
invaluable to test the validity of the models, if they were made 
in a manner that allov;ed us to get simultaneous rainfall rate 
data. If the problem of getting enough aircraft data to compare 
with the SMMR satellite data is severe {primarily, or course, 
this is due to financial considerations) , the problem of getting 
rain rate data on the ground (or rather the sea surface) from 
rains gagues is worse. There is no network of rainfall gagues 
on the ocean, partly for economic reasons, and partly becasue 
it is difficult to make the measurement with existing types of 
equipment while at sea. The only available data would have to 
come from the few radar weather stations that are close enough 
to the shore to provide adequate coverage over the ocean, e.g., 
the one in Miami, Florida. It might also be possible to use 
ship-bourn radar, if such a facility could be made available. 
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APPENDIX A - ATMOSPHERE OPACITY 


In Table A-I, we show the amount of water vapor and 
the no-rain optional depth at each frequency for the models 
used in this study. 
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TABLE A-1 



ATMOSPHERIC OPACITY TABLE FOR A LAPSE RATE OF 6.5 deg/km 




a t m o 

s p h e r 

i c op 

a c i t y 


Freezing 

Water 

.22 

.356 

.6 

.7 

1.23 

Level 

(km) 

Vapor 

(cm) 

-1 

cm 

-1 

cm 

-1 

cm 

-1 

cm 

-1 

cm 

5 

7.056 

1.3679 

-2 

X 10 

2.8707 

-2 

xlO 

.1365 

.3686 

.2090 



1.2235 

2.3682 

9.8791 

.2504 

.1744 

4 

4.534 

X lO"^ 

X 10“^ 

xlO”^ 





1.13 

2.0196 

7.2203 

.1666 

.1495 

3 

2.715 

-2 

X 10 

-2 

X 10 

-2 

X 10 





1.0736 

1.782 

5.36 

.1088 

.1318 

2 

1.567 

-2 

X 10 

-2 

X 10 

-2 

X 10 
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